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In models of low-energy gauge mediation, the observed Higgs mass is in tension with the cosmo¬ 
logical limit on the gravitino mass 7113/2 16 eV. We present an alternative mediation mechanism 

of supersymmetry breaking via a 17(1) D-term with an Ae-inspired particle content, which we call 
vector mediation. The gravitino mass can be in the eV range. The sfermion masses are at the 10 TeV 
scale, while gauginos around a TeV. This mechanism also greatly ameliorates the /r-problem. 


Introduction. Supersymmetry (SUSY) is the prime 
candidate for solving the naturalness problem of the 
Standard Model Higgs sector (see, e.g., [I])- Models of 
low-energy SUSY breaking are particularly attractive be¬ 
cause the light gravitino G can also be produced at col¬ 
liders with the striking signature of additional photons in 
the decay of the bino B -T Gy [2] • The most extensively 
studied mechanism for this is gauge mediation [3]. 

However, cosmological limits on the gravitino are in 
strong tension with low-energy gauge mediation. In the 
Minimal Supesymmetric Standard Model (MSSM), the 
observed Higgs mass of 125 GeV HIS] requires large A- 
terms or stop masses in the 10 TeV range. In gauge me¬ 
diation, generating large A-terms through renormaliza¬ 
tion group evolution [ 6 ] or large stop masses requires the 
SUSY breaking scale to be above 1000 TeV, and hence 
the gravitino mass TO 3/2 to be above a keV [7]. Such a 
gravitino would overdose the Universe unless the reheat¬ 
ing temperature after the inflation is kept unnaturally 
low 019], leaving little room for any conceivable baryo- 
genesis scenario. When TO 3/2 < 1 keV, the gravitino 
does not overdose the Universe, but would be a warm 
dark matter component, ruled out by the Lyman-a for¬ 
est data unless TO 3/2 ;< 16 eV [TU]. The tension may be 
eased if the Higgs sector is extended to the NMSSM or 
Dirac-NMSSM if the messenger sector couples to 
the Higgs [12], or is strongly coupled However, 

these models are somewhat more complicated. 

In this Letter, we present an alternative simple medi¬ 
ation mechanism of SUSY breaking, where the sfermion 
masses are generated at tree-level from a t/(l) H-term. 
TO 3/2 can be brought down to the eV scale, solving the 
cosmological problem. In addition, the gluino is most 
likely within the reach of the HL-LHC. 

Vector Mediation We first present the basic idea of 
what we call the vector mediation of SUSY breaking. 

We assume that the SUSY breaking sector produces 
both an F-term and a F-term of a [/(I) vector multi- 
plet under which the Standard Model matter fields are 
charged. In order for all of the sfermions in the Standard 
Model to acquire positive soft mass-squared, all quark 


and lepton multiplets must have U{1) charges of the 
same sign (positive). Anomaly cancellation requires ad¬ 
ditional superfields with negative U(l) charges, which ac¬ 
quire negative soft mass-squared from the F-term. They 
must come in Standard Model vector-like representations 
so that the vector-like masses can overcome the negative 
soft mass-squared. It is economical if the vector-like su- 
permultiplets also act as messengers that generate gaug- 
ino masses at the one-loop level. Because of the tree-level 
generation of scalar masses, the scalar masses as well as 
the scale of SUSY breaking are both in the 10-100 TeV 
range, solving the cosmological gravitino problem. 

The tree-level mediation of SUSY breaking via a F(l) 
F-term was also considered in Refs. mnnj. Their inter¬ 
est was in a high-scale mediation mechanism which has a 
much more weakly coupled goldstino and thus have very 
different cosmology and collider signatures [l9] . 

An Explicit Model The simplest particle content that 
satisfies the requirements explained above is one inspired 
by Eq unification. By embedding S'O(IO) x F(l),/, into 
Eq, the fundamental representation 27 decomposes as 

27 = T'(16, -f 1) 0 $(10, -2) 0 5(1, 04). (1) 

We identify $ as the quark and lepton supermultiplets, 
while $ plays the role of messengers as well as the Higgs. 
S is responsible for the SUSY breaking. 

SUSY Breaking We introduce a vector-like mutiplet 
A(—4) and Y(04) charged under U{1)^. We also intro¬ 
duce a neutral field Z{0) and write the superpotential 

W = MSX + XZ{XY -v^). (2) 

It is essential that the model here is chiral under F(I)^, 
since charge conjugation invariance would prevent the 
generation of a F-term. Note that the above superpo¬ 
tential is generic and F(l)/j invariant, which guarantees 
that it breaks SUSY |2D| [31] . 

This model breaks SUSY with a positive F-term, 
which induces tree-level positive mass-squared to the 
matter fields $. Motivated by models of dynami¬ 
cal SUSY breaking, a la Izawa-Yanagida-Intriligator- 
Thomas (lYIT) [221 US], we take the values A = 47 r, 
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e = 0.1 and M = 2v as an example. We will use the 
notation Xq = {X), etc. The minimum is at Xq = 0.58u, 
Yq = 1.71u, S'o = 0 and Zq = 0. The non-zero SUSY 
breaking parameters are Fs = MXq = 1.16u^, Fz = 
X{XoYo - v^) = -0.1^2, and D = 4e(-|Yo|2 + |YoP) = 
l.02v^ > 0. All particles here are massive except for the 
goldstino, which is a linear combination of the S fermion 
and the U(l)^ gaugino. Hence there is no cosmological 
Polonyi problem |24) . 

Note that the size of the H-term is completely fixed by 
the size of the F-term as 

D=^F^{4e)Fs, = 2(4e)2(|Xop + lYoH. (3) 

Here, My is the mass of the U{1)^ gauge boson. This 
equation is a simple consequence of the equation of mo¬ 
tion for the massive vector multiplet. 

In our Eq inspired model, there are two additional S' 1.2 
fields which do not play a role in SUSY breaking. To give 
S' 1.2 a mass, we introduce two new neutral fields N 12 
and introduce the couplings W = kA^i. 2 Si. 2 Y. Once X 
obtains its vacuum expectation value, this superpotential 
combined with the F-term gives a mass to Si .2 and Ni^ 2 - 

Tunneling to the true vacuum In order to generate 
gaugino masses at one-loop, we need the superpotential 

IU„ =^S$$. (4) 


Because this superpotential does not respect the U(l)u 
symmetry, a supersymmetric vacuum appears. It is lo¬ 
cated at 


X - ^ 
Z-o- y 

ro 


^^0 z ^ 9 M ^ _ 2gMv'^ 
2M ’ ° fcA ’ ° 


(5) 


where the value of $0 is determined by minimizing the 
F-term. Because there is both a SUSY-preserving and 
a SUSY-breaking vacuum, one must consider tunneling 
between the two vacua. From the form of the supersym¬ 
metric vacuum, it is clear that if k is small, the distance 
between the two vacua will be large. In the limit that 
k ^ 0, the SUSY vacuum disappears and the metastable 
vacuum becomes stable. 

The life-time of the metastable vacuum can be cal¬ 
culated using semiclassical techniques |25j . The tun¬ 
neling rate per unit volume is T oc e~^ where B = 
Se{ 4 >{i")) — ^EiM- 4 ’o is the field configuration of the 
metastable vacuum, Se is the Euclidean action and 4){r) 
is the bounce profile. The Euclidean action is 


SE{(t>{r)) = 27r2 





v{m) , 
( 6 ) 


where the sum goes over all five fields, (j)i = 
{X, Y, Z, 4), S). (j){r) is the solution to 


d(j)i dV 

dj-'i j- dj- d(l)i ’ 


Bounce action 



FIG. 1. Contours of the vacuum bounce action on the {k,g) 
plane. The solid blue, red, black, purple and gray lines are 
for B = 250, 350, 450, 550, and 650, respectively. For the 
universe to be stable, B > 450. The dashed black line and 
the grey region show where the SUSY breaking minimum be¬ 
comes tachyonic. We expect that the approximations used to 
calculate the bounce action fail in this region. 


subject to the boundary conditions 

= 0) =0, 0i(r^oo) = ^o- (8) 

We need B > 450 for the metastable minimum to live 
longer than the age of the universe. 

Calculating the bounce properly requires solving a five¬ 
dimensional differential equation. To simplify the com¬ 
putation, we approximate the tunneling as confined to a 
one-dimensional sub-space proceeding in a straight line 
between the two minima. In this approach, solving for 
the bounce action is numerically very stable. 

We expect this approximation to be valid when the 
potential has large first derivatives, e.g., when the maxi¬ 
mum separating the two minima is very large. Then the 
friction term, the I/r term in Eq. 0. is important. In 
order for the bounce to travel between the two minimum, 
it first stays near the SUSY minima until large r where 
the friction term is irrelevant. The bounce then rolls to 
the SUSY breaking minimum. In this limit, because the 
bounce action is required to start very close to the SUSY 
minima and to end at the SUSY breaking minimum, we 
expect that this one-dimensional approximation is good. 
We have also checked that the maximum separating the 
two minima lies on the straight line between the two min¬ 
ima to within a few percent. 

In Eig. we plot contours of the bounce action on the 
{k, g) plane. The shaded region is where the messengers 
go tachyonic. Then there is no maximum between the 
two vacua and the bounce does not start near the SUSY 
vacua. Thus, we expect the approximation used to calcu¬ 
late the bounce action to be very poor. We have verified 
numerically that the one-dimensional approximations of 
the bounces in these regions of parameter space do not 
start near the SUSY vacua. 

Aside from the vacuum transition rate, there is also the 
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FIG. 2. Contours of on the {k,g) plane. The solid 

blue, red, and black lines are for = 130, 230, and 330, 

respectively. For the universe to be stable, > 230. 

The high temperature expansion is used to estimate Tc and 
hence the result is meant to be only qualitatively correct. See 
also the caption of Fig. 


probability of thermally tunneling between the two vacua 
in the early universe. For T > v, the SUSY vacuum dis¬ 
appears and there is only the SUSY-breaking vacuum. 
At a lower temperature, the SUSY-preserving vacuum 
appears. At the temperature Tc, the SUSY vacuum be¬ 
comes the true minimum and tunneling between the two 
vacua can occur. We obtain the thermal bounce action, 
Bth = S'th((/>(?’)) - -S'th(</>o), by solving the three dimen¬ 
sional versions of Eqs. ® ■ 

To evaluate the thermal tunneling rate F ^ 
we make the following approximations. We first assume 
that the thermal bounce Bth is a constant and deter¬ 
mine it numerically using the zero temperature potential, 
again using the one-dimensional approximation. We find 
that in order for the universe to be stable, Bth/Tk > 230. 
The second approximation we use is the high temperature 
expansion to determine T^. We make this approximation 
because it is numerically difficult to find when minima 
appear in a five dimensional space when the exact one 
loop integrals do not converge quickly. We find Tc tends 
to be smaller than v, where the high temperature expan¬ 
sion is not valid. Therefore, the result here should be 
viewed as only qualitatively correct. 

We plot the contours of constant Bth/Ii, on the {k,g) 
plane in Fig. If the value of Bth/Tc were off by a 
factor of two, the limit on k changes by 20 % for a fixed 
g. Again, the one-dimensional approximation is poor in 
the shaded region where messengers go tachyonic. 

SUSY Spectrum All of the squarks and sleptons ob¬ 
tain the same mass-squared from the non-zero B-term, 


m|, = eB > 0. (9) 

and hence there is no flavor problem. 

As mentioned earlier, gaugino masses are generated 
from one-loop diagrams involving <&. The mass of the 
fermion component of <i> is = gYo, while the boson 


components have the mass matrix (see Eq. 


1 


£3 


f (gYof - 2eB kFg 
I kFs ( 5 Yo )2 - 2eB 


$ 

( 10 ) 

We choose g = 0.5 and k = 0.35, in addition to the pa¬ 
rameters in the previous section. There are no tachyons 
and the metastable vacuum is stable on timescales of or¬ 
der the age of the universe. 

At this point we need to discuss the Higgs. In super- 
symmetric Eq theories, we can regard the doublets in $ as 
the Higgs superfields. We assign positive matter parity to 
doublets in one of $ so that it can have Yukawa couplings. 
On the other hand, all triplets and other doublets have 
negative matter parity so that they do not have Yukawa 
couplings to Ik. They are simply messengers for the gaug¬ 
ino mass. The Higgses have the same mass matrix as in 
Eq. (10 1 , and we fine-tune p, = gnYo against 2eB to ob¬ 
tain a mass-squared at the 100 GeV scale. It represents 
a tuning at the level of 10“^. Bp = knFs can also be 
made at the 100 GeV scale by choosing k^ « 10“"^. Note 
that the bare /i-parameter is forbidden and is generated 
only with B(l),/, breaking, ameliorating the /i-problem 
from the Planck scale to the 100 TeV scale. 

We have three color triplets and two electroweak 
doublets acting as messengers to generate the gaugino 
masses. The 1-loop results for the gaugino masses can 
be found in Eq. (2.3) of Ref. [5^. The gaugino masses 
can be enhanced relative to the standard gauge media¬ 
tion 3^ due to the B-term. Eor the above numerical 
example, the enhancement factor is 1.39. The on-shell 
gluino mass is further enhanced from M 3 by the QGD 
correction by 1 + | _ 

In these models, the scalars are heavy, the gauginos 
are light and the gravitino is the Lightest Supersymmet¬ 
ric Particle. Eor this type of mass spectrum, the current 
ATLAS bound on the gluino mass is 1.28 TeV [27] while 
the current CMS bound is around 1.15 TeV [28]. Because 
our gaugino masses are enhanced with respect to stan¬ 
dard gauge mediation, it is not very difficult to obtain 
a heavy enough gaugino. Eor the choice of parameters 
mentioned earlier, we find v > 87 TeV from the gluino 
mass limit. 

In the large A limit, Y and Z can be integrated out, 
and the mass spectrum of the singlets can be computed 
as follows. Among the singlets, S' is a complex scalar and 
X has only its real part (the imaginary part is eaten by 
the U{1)^ gauge boson), with masses 


TOg = 


XI 


x! 


Yi 


+ l(ye\-XlYYS), ( 11 ) 


m 


X 


= M 


2A2 + 16e2(3A2 + 5Yo2) 


^ 0 ^ 




( 12 ) 


Eermions of X and S become a Dirac particle of mass 


m% = M" 


XS 


X^ + Yi 


+ 32e^(X^ + Yq^). (13) 
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FIG. 3. An example mass spectrum in our model. The fields 
are presented according to the S'O(IO) representations, 1, 10, 
and 16, as well as the Standard Model gauginos. The grav- 
itino is too light to be shown. 


The C/(l)^ gauge boson acquires a mass 

m^ = 32e2(x2 + Fo2), 


(14) 


We have already presented the mass spectrum of the 4) 
and 4/ multiplets. We verified that StrM^ = 0 once spec¬ 
trum integrated out at 0 (A^(Xq -|- Yq )) is included. 

In order for the gluino to be heavier than a TeV, we 
typically have scalars at 20-100 TeV. If e ~ 0(1), this 
is a similar mass spectrum to mini-split |29j . pure grav¬ 
ity mediation [30] . and SUSY breaking mechanisms with 
accidental ii-symmetries, e.g.^ metastable Intriligator- 
Seiberg-Shih (ISS) SUSY breaking [3T]. The main dif¬ 
ference of vector mediation with previous models is the 
presence of a superlight gravitino. The mass spectrum 
is shown in Fig. The radiative corrections from stops 
at the 10 TeV scale easily raise the lightest Higgs boson 
mass to the observed mass of 125 GeV. Note that the fine- 
tuning is smaller compared to other models with a split 
spectrum because the mass of three 16 is at ~ 30 TeV 
whereas the other scalar masses are at ~ 100 TeV. This is 
because the U{1).^ charges differ by a factor of four and 
the H-term is always smaller than the F-term. These 
scalars may be discovered at a 100 TeV pp collider. 

Implications The most important implication is that 
the gravitino mass is 


m 


2 

3/2 


V 


3Mpi 


1.9v^ 

iMpi 


(2.5 eV)2 


(15) 


for the above parameters. Here, Mpi = 2.4 x 10^® GeV. 
Such a low mass gravitino is not possible in the usual 
low-energy gauge mediation and is completely compatible 
with the cosmological limit of 77 x 3/2 < 16 eV. There is a 
clear upper limit on the gluino mass from the cosmology 
of about 5 TeV. The LHG would most likely discover 
the gluino, and a 100 TeV pp collider would completely 
close the window. The bino decays promptly B —> 7 G, 


giving evidence of low-energy SUSY breaking. The ratio 
among gaugino masses is also unique, different from the 
standard gauge-mediation or unification scenarios. 

The right-handed neutrino is charged under U (1).^ and 
hence we cannot use the seesaw mechanism for neutrino 
masses. One can use a linear combination of U{\)^ and 
U{l)y^ instead to make the right-handed neutrino neutral 
[32] . This is an interesting variation to be considered. In 
the simplest realization of these models presented in this 
Letter, the right-handed neutrinos are Dirac, are present 
at low energies, and can contribute to energy densities 
at BBN. Much of the number density in right handed 
neutrinos is diluted away by the entropy generated at 
the QCD phase transition, resulting in ~ 0.19. This 
may be verified in the future cosmological data. 
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